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Fast thermalization and a strong buildup of elliptic flow of QCD matter as found at RHIC are 
understood as the consequence of perturbative QCD (pQCD) interactions within the 3+1 di- 
mensional parton cascade BAMPS. The main contributions stem from pQCD bremsstrahlung 
2^3 processes. By comparing to Au+Au data of the flow parameter 112 the shear viscosity 
to entropy ratio rj/s has been extracted dynamically and lies in the range of 0.08 and 0.2. Also 
jet-quenching has been investigated consistently within a full dynamical picture of the heavy 
ion collision. The results for gluonic jets indicate a slightly too large suppression, but are 
encouraging to understand the two major phenomena, strong flow and jet-quenching, within 
a unified microscopic treatment of kinetic processes. In addition, simulations on the temporal 
propagation of dissipative shock waves lead to the observation that an ri/s ratio larger than 0.2 
prevents the development of well-defined shock waves on timescales typical for ultrarelativistic 
heavy-ion collisions. 

1 Overview 

The large elliptic flow parameter V2 measured by the experiments at the Relativistic Heavy Ion 
Collider (RHIC )^ suggests that in the evolving QCD fireball a fast local equilibration of quarks 
and gluons occurs at a very short time scale < 1 fm/c, and that the locally thermalized state of 
matter created, the quark gluon plasma (QGP), behaves as a nearly perfect fluid. In addition, 
the phenomenon of jet-quenching has been another important discovery at RHIC^. So far, it 
has not been possible to relate both phenomena by a common understanding of the underlying 
microscopic processes. 

In this talk, we demonstrate that perturbative QCD (pQCD) collisional and, more impor- 
tantly, radiative interactions can explain a fast thermalization of the initially nonthermal gluon 
system, the large collective effects of QGP created at RHIC and the smallness of the shear 
viscosity to entropy ratio in a consistent manner by using a relativistic pQCD-base d on-shell 
parton cascade Boltzmann approach of multiparton scatterings 
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Also, due 

to the inclusion of radiati ve p rocesses. BAMPS can simultaneously account for the quenching 
of high- momentum gluons In addition, the possible propagation of dissipative shock waves 
in the QGP can be studied^!!'. 

2 Elliptic Flow and Shear Viscosity 

BAMPS is a parton cascade, which solves the Boltzmann equation and is based on the stochastic 
interpretation of the transition rate"^, which ensures full detailed balance for multiple scatterings. 
Gluon interactions included are elastic and screened Rutherford-like pQCD gg gg scatterings 




Figure 1: (Color online) Left panel: Elliptic flow vs. the number of participating nucleons for Au+Au collisions 
at y^SNN ~ 200 GeV. The points are STAR and PHOBOS data for charged hadrons within the rapidity |y| < 0.5 
and \y\ < 1, respectively, whereas the curves with symbols are results for gluons within \y\ < 1, obtained from the 
BAMPS calculations with Ua = 0.3 and 0.6 and with two freezeout energy densities, ec = 0.6 and 1 GeVfm"^. 
Right panel: the shear viscosity to entropy density ratio rj/s at the central region during the entire expansion. 
rj/s values are extracted from the simulations at impact parameter b = 4, 6.3, and 8.6 fm. The upper band shows 
the results with Os = 0.3 and the lower band the results with Us = 0.6. 



as well as pQCD inspired bremsstrahlung gg ^ ggg of Gunion-Bertsch typeP^^^^^. In the default 
simulations, the initial gluon distributions are taken in a Glauber geometry as an ensemble of 
minijets with transverse momenta greater than 1.4 GeV^, produced via semihard nucleon- 
nucleon collisions. The minijet initial conditions and the subsequent evolution using the present 
prescription of BAMPS for two sets of the coupling a<j = 0.3 and 0.6 give nice agreements 
to the measured transverse energy per rapidity over all rapidities'^^. Other initia l con ditions 
incorporating gluon saturation such like the color glass condensate are also possible^^El^ 

The elliptic flow V2 = {{Px —Py)/PT) be directly calculated from microscopic simulations 
fSEl. The results are compared with the experimental data, assuming parton-hadron duality. 
The left panel of Fig. 1 shows the elliptic flow V2 at midrapidity for various centralities (impact 
parameters). Except for the central centrality region either the results with as = 0.3 and Cc = 0.6 
GeVfm"^ (green curves with open triangles) or those with as = 0.6 and Cc = 1 GeVfm"'^ (red 
curves with open squares) agree perfectly with the experimental data^^^. Within the Navier- 
Stokes approximation the shear viscosity ry is directly related to the so called transport rate 
^ and can be calculated dynamically and locally in a full and microscopical simulation (see 
the right panel of Fig. 1). Bremsstrahlung and its back reaction lower the shear viscosity to 
entropy density ratio rj/s significantly by a factor of 7, compared with the ratio when only elastic 
collisions are considered. According to the extraction by means of the full simulati on giv en in 
Fig. 1, rj/s is most probably lying between 0.2 for ag = 0.3 and 0.08 for as = 0.6^^^. The 
latter value matches the lower bound of tj/s = \/Att from the AdS/CFT conjecture. The shear 
viscosity slightly increases up to 20% when calculating it by means of Grads momentum method 
and using the 2nd order Israel-Stewart framework of dissipative relativistic hydrodynamics^^. 



3 Jet-Quenching and Dissipative Shock Waves 

Jet-quenching is generally specified in terms of the nuclear modification factor Raa- We directly 
compute Raa by taking the ratio of the final pT spectra to the initial mini-jet spectra. For this 
a suitable weighting and reconstruction scheme of initial jet-like partons had to be developed 
Figure 2 shows the result for the gluonic contribution to Raa-, exhibiting a clear suppression of 
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Figure 2: (Color online) Gluonic Raa at midrapidity (i/ e [—0.5, 0.5]) as extracted from simulations for central 
Au+Au collisions at 200 AGeV. The shaded area indicates the statistical error. For direct comparison the result 
from Wicks et al. for the gluonic contribution to Raa and experimental results from PHENIX for and STAR 

for charged hadrons are shown. 

high-pT gluon jets at a roughly constant level of R^^^^ ^ 0.053. 

The suppression of gluon jets is approximately a factor 3-^4 stronger than seen in exper- 
imental pion data^Sl. An excessive quenching, however, was to be expected since at present 
the simulation does not include quarks, which are bound to lose considerably less energy due to 
their color factor and dominate the initially produced jets from px ~ 20 GeV. Indeed, comparing 
with state of the art results from Wicks et al. for the gluonic contribution to Raa (seen as 
the line in Fig. [2]), which in their approach together with the quark contribution reproduces 
the experimental data, one finds better agreement. The inclusion of light quarks will provide 
important means of further verification and will be addressed in a future work. 

In the context of jet-quenching, exciting jet-associated particle correlations have been 
observed, which indicates the formation of shock waves in form of Mach cones induced by 
supersonic partons moving through the quark-gluon plasma. Shock waves can only form and 
propagate if matter behaves like a fluid. At present it is not known whether the ij/s value 
deduced from V2 data is sufficiently small to allow for the formation of shock waves. A step 
towards this problem is to consider the relativistic Riemann problem in viscous gluon matter^. 

In Fig. 3 we show the solution of the Riemann problem for various rj/s values as computed 




Figure 3: (Color online) The pressure (left) and the velocity (right) profiles at time t = 3.2 fm/c. At t = 0, the 
pressure is Pq = 5.43 GeVfm"^ for z < and P4 = 0.33 GeVfm"^ for z > 0. 



with BAMPS. The results demonstrate the gradual transition from the ideal hydro dynamical 
limit to free streaming of particles. A larger r]/ s value results in a finite transition layer where 
the quantities change smoothly rather than discontinuously as in the case of a perfect fluid. 
A nonzero viscosity smears the pressure and velocity profiles and impedes a clean separation 
of the shock front from the rarefaction fan. A criterion for a clear separation, and thus the 
observability of a shock wave, is the formation of a shock plateau, as in the ideal fluid case. The 
formation of a shock plateau takes a certain amount of time. From Fig. 3 we infer that, for 
rj/s > 0.1, a shock plateau has not yet developed at t = 3.2 fm/c, whereas for rj/s < 0.1, it has 
already fully formed. Applying time scaling arguments, the formation of shock waves in gluon 
matter with ij/s > 0.2 probably takes longer than the lifetime of the QGP at RHIC — . 

In summary, within the pQCD based parton cascade BAMPS early thermalization, the build 
up of a large elliptic flow V2, the smallness of the rj/s ratio, and the quenching of gluonic jets in 
Au+Au collisions at RHIC can be understood within one common setup. This is a committed 
and large scale undertaking. Further analyses on jet quenching, on particle correlations, on 
possible Mach cone like behaviour initiated by a jet, on quark degrees of freedom including 
hadronization, on initial conditions and on the exploration and use of dissipative hydrodynamics 
are underway to establish an even more global picture of heavy ion collisions. 

Acknowledgments 

This work was supported by BMBF, DAAD, DFG, GSI and by the Helmholtz International Cen- 
ter for FAIR within the framework of the LOEWE program (Landes-Offensive zur Entwicklung 
Wissenschaftlich-okonomischer Exzellenz) launched by the State of Hesse. 

References 

1. S. Adler et al. (PHENIX Collaboration), Phys. Rev. Lett. 91 182301 (2003); J. Adams et 
al. (STAR Collaboration), ibid. 92 052302 (2004). 

2. C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 89, 202301 (2002); K. Adcox et 
al. (PHENIX Collaboration), ibid. 88, 022301 (2002). 

3. Z. Xu and C. Greiner, Phys. Rev. C 71 064901 (2005). 

4. Z. Xu and C. Greiner, Phys. Rev. C 76 024911 (2007). 

5. A. El, Z. Xu , and C. Greiner, Nucl. Phys. A806, 287 (2008). 

6. Z. Xu, L. Cheng et al., larXiv:0812.3"839] [hep-ph]. 

7. Z. Xu and C. Greiner, Phys. Rev. Lett. 100, 172301 (2008). 

8. Z. Xu, C. Greiner, and H. Stocker, Phys. Rev. Lett. 101, 082302 (2008). 

9. Z. Xu and C. Greiner, Phys. Rev. C 79, 014904 (2009). 

10. O. Fochler, Z. Xu , and C. Greiner, Phys. Rev. Lett. 102, 202301 (2009). 

11. I. Bouras et al., arXiv:0901T927j [hep-ph]. 

12. J. Adams et al. (STAR Collaboration), Phys. Rev. C 72 014904 (2005); B. B. Back et al. 
(PHOBOS Collaboration), ibid. 051901 (2005). 

13. A. El et al., Phys. Rev. C 79 044914 (2009). 

14. J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 91, 172302 (2003); A. Adare et 
al. (PHENIX Collaboration), ibid. 101 232301 (2008). 

15. S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, Nucl. Phys. A784, 426 (2007). 

16. F. Wang (STAR Collaboration), J. Phys. G 30, S1299 (2004); J. Adams et al. (STAR 
Cohaboration), Phys. Rev. Lett. 95, 152301 (2005); S. S. Adler et al. (PHENIX Collab- 
oration), ibid. 97, 052301 (2006); J. G. Ulery (STAR Collaboration), Nucl. Phys. A774, 
581 (2006); N. N. Ajitanand (PHENIX Collaboration), ibid. 783, 519 (2007); A. Adare et 
al. (PHENIX Collaboration), Phys. Rev. C 78, 014901 (2008). 



